Single-qubit gates in two steps with rotation axes in a single plane 
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Any single-qubit unitary operation or quantum gate can be considered a rotation. Typical ex- 
perimental implementations of single-qubit gates involve two or three fixed rotation axes, and up to 
three rotation steps. Here we show that, if the rotation axes can be tuned arbitrarily in a fixed plane, 
then two rotation steps are sufficient for implementing a single-qubit gate, and one rotation step 
is sufficient for implementing a state transformation. The results are relevant for "exchange-only" 
logical qubits encoded in three-spin blocks, which are important for universal quantum computation 
in decoherence free subsystems and subspaces. 
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I. INTRODUCTION 

In the quantum circuit model [I] , a universal quantum 
computer requires an entangling two-qubit gate such as 
CNOT, and a set of single-qubit gates [2 j. Although a 
finite set of quantum gates is sufficient for universality, 
fault-tolerant circuits require a large number of gates, 
even for very simple operations [3 . It is therefore impor- 
tant to be able to perform gate operations as efficiently as 
possible, in order to minimize the effects of decoherence 
or gating errors. 

Single-qubit operations can be viewed as rotations of 
the state vector of a logical qubit on a unit Bloch sphere 
[3]. The most efficient method for rotating a spin qubit 
would be to apply a magnetic field along the desired axis 
of rotation; however this is not practical for most physi- 
cal implementations. For example, to independently con- 
trol an array of electron spin qubits in quantum dots [4 
would require an array of tunable micromagnets, which 
is experimentally challenging. 

A more common approach is to provide two or three 
fixed, orthogonal rotation axes. This enables arbitrary 
rotations in three or fewer steps, for example, by using 
the Euler angle construction. Many qubit implementa- 
tions employ this strategy. For single- spin qubits, this 
could involve tunable, external magnetic fields oriented 
along two orthogonal axes, although this could still be 
challenging in a scalable architecture. For logical qubits 
composed of two or more physical spins, it is possible 
for different physical coupling mechanisms to control the 
different rotation axes on the Bloch sphere. For exam- 
ple, singlet-triplet logical qubits formed in double quan- 
tum dots use local magnetic field gradients to generate 
rotations about one axis, and exchange interactions to 
generate rotations about an orthogonal axis [5]. 

In this paper, we are concerned with intermediate situ- 
ations, between these two extremes, where rotations can 
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FIG. 1. (Color online) Single-plane rotation method: the 
shaded regions indicate the xz plane of the Bloch sphere, 
and dashed blue arrows indicate the rotation axes, (a) An 
arbitrary transformation, from the state \^a) along n^, to 
the state \iPb) along n^, via a single rotation about the axis 
n. In this example, nA = (0, 1, 1)/V% and n.B=(l, 1, 0)/\/2. 
As described in the main text, we obtain the rotation 
axis n=(l, 0, l)/v / 2, and rotation angle <f> — —cos -1 (1/3). 
(c) An arbitrary single qubit rotation, performed in two 
steps. In this example, 7Z(y, 37r/2)=7^(n2, ?r)7^(z, tt) with 
n 2 =(-l,0,l)/V2. 



be performed about an arbitrary axis constrained to lie in 
a single plane. We consider the example of an exchange- 
only logical qubit in a triple quantum dot [7] , which is 
the main motivation for our work. Exchange-only qubits 
are of interest because they are formed in decoherence 
free subspaces and subsystems [8UTU]. and they can be 
operated using only fast exchange interactions [TTHT3] . 
A similar, single-plane rotation scheme is also possible 
for chirality-based logical qubits in a triple dot [14]. In 
fact, for any system where the effective magnetic fields 
along two orthogonal axes can be tuned arbitrarily and 
simultaneously, the single-plane rotation condition is sat- 
isfied. 
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We will build upon the known result |6 that the inter- 
nal couplings between the physical spins in an exchange- 
only qubit can be used to generate a continuous set of ro- 
tations in the xz plane of the logical qubit. We will show 
that such single-plane rotations reduce the total number 
of steps required for single-qubit operations. For exam- 
ple, a state transformation, which maps a specific ini- 
tial state onto a specific final state, can be accomplished 
in one step, while an arbitrary single-qubit rotation can 
be accomplished in two steps. We provide constructive 
proofs for both of these problems. 

II. SINGLE-STEP STATE TRANSFORMATION 

We first prove that a given initial state |^) can be 
transformed to a specified state \^b), by a single rota- 
tion about an axis in the xz plane. The Bloch sphere 
geometry is shown in Fig. [TJa). Specifically, we want to 
determine the rotation axis n and the rotation angle <p 
that satisfy 

ll(n,<p)\* A ) = \* B ) , (1) 

where the rotation operator 7£(n, 0)=exp(— icr • n</>/2) is 
defined in terms of the Pauli spin matrix vector a. The 
qubit state vectors and are also pictured on the 
Bloch sphere of Fig. [TJa). If the rotation axis is allowed 
to point in any direction (not just the xz plane), then we 
could choose n oc (n^ + n#) with <j) = tt. We now show 
that the desired result can be achieved, even when n is 
confined to the xz plane. 

To perform a state transformation in a single step, it 
is clear that the rotation axis must be equidistant from 
both ua and n^. This constraint defines a plane, given 
by n • (ua - n^) =0. On the other hand, we require the 
rotation axis to lie in the xz plane, which is defined by 
n • y = 0. The intersection of the two planes is given by 

n= V X ^- AB \ r (2) 
|y x (n A - n B )\ 

Figure [TJb) shows the projection of and in the 
plane perpendicular to n. The inscribed angle is the angle 
of rotation, 0. Since the length of all three sides of the 
triangle are known, as indicated in the figure, we can 
obtain 

cos(/>=l-|n A -n B | 2 /2L 2 , (3) 

where L = \ua — (n • £u)n| = \h B — (n • ri£)ri|. The sign 
of (j> is given by sgn[(n^ x n#) • n]. When the two planes 
are equivalent, we cannot use Eq. However in this 
special case, + lies in the xz plane, and we can 
simply choose n ex (n^ + rig) with (j) = tt. 

Figure [TJa) shows an example of the state trans- 
formation procedure for the case n^=(0, ^) and 
n B =(^,^,0). From Eqs. ^ and we obtain 
the rotation axis n=(^,0, ^) and the rotation angle 
-cos-^l/S) ~ -70.53°. 
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FIG. 2. (Color online) Definition of the angular variables, 
(a) The logical rotation axis n, with the polar angle and 
the azimuthal angle ip. (b) The projection of n onto the xz 
plane. We define a new coordinate system with axes x' and 
z', such that n' x — 0, by rotating the xz plane around the y 
axis by angle /3. 

III. TWO-STEP QUBIT ROTATIONS 

We now provide a constructive proof that any single- 
qubit gate (up to a global phase) can be generated in 
two rotation steps, when the rotation axes point in an 
arbitrary direction in a single plane. Specifically, we want 
to solve for the rotation parameters defined by 

K(h, <P) = e^TZ 2 (n 2 , 2 )ftx (ni ,4>i) , (4) 

where the rotation axis n can point anywhere in the 
Bloch sphere, but the individual rotation axes ni and 
n2 lie in the xz plane. It is convenient to work with 
angular coordinates defined by 

n = (sin # cos?/;, sin # sin cos#) , (5) 
ni = (sin#i, 0, cos#i) , (6) 
n 2 = (sin #2,0, cos # 2 ) , (7) 

and illustrated in Fig. |2|a). Since an arbitrary ro- 
tation is characterized by three parameters (0, -0, 0), 
while the right- hand- side of Eq. Q involves five pa- 
rameters (??, #1, (/>i, #2, ^2), the equation is clearly under- 
constrained; many solutions exist, any of which suits our 
needs. 

We now demonstrate that at least one solution exists 
by providing an explicit, analytical construction. We first 
simplify the problem by transforming to a new set of 
coordinate axes defined by the projection of the logical 
rotation axis n onto the xz plane. The normalized pro- 
jection becomes our new z axis, as shown in Fig. |2|b). 
The rotation angle for the transformation, (3, is given by 

cos /3 = Uz — , sin /3 = Ux — . (8) 

V n l+ n2 z V n<2 x+ n2 z 

where n x and n z are the components of n in the original 
coordinate system. Using primed variables to indicate 
the new coordinate system, we have 

n = sinfl'y + cos(9 / z / . (9) 
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where 



sin 0' = sin sin ip . 



(10) 



Note that y is the same in both coordinate systems. 

We now solve Eq. Q in the primed coordinate sys- 
tem. Without loss of generality, we will restrict the two 
rotation axes in the range 0[, 6' 2 G [0, tt) and the rotation 
angles in the range 0i, 02 G [0, 27r), since a rotation with 
0\ > tt by an angle & is equivalent to the rotation with 
• — tt by an angle 27r — ^ . We can take advantage of the 
under-constrained nature of the problem by making the 
convenient choice </>2 = tt. (See Appendix [B] for details.) 
We then find 







k cos sin - 



k sin # tan - 



cos - 



tan 0[ , 



(ii) 

(12) 
(13) 

and 
re- 
spectively. These quantities are related to the original 
coordinate system through 6i=6[ + /? and #2 =#2 + 

As noted above, Eqs. ( 11 )-( 13) do not represent unique 
solutions to Eq. Q. For example, we may obtain an 
alternative solution with the choice <p\ = tt. After a 
similar analysis, we then obtain 



where k = sgn^] = sgn^in^sin^]. Here 
0[ are determined from Eqs. (11), (p|), and (13) 



y 2> 



e\ =0 



-k cos sin - 



fc sin 6' tan 



cos 



tan #9 . 



(14) 
(15) 

(16) 



To provide a practical comparison, we consider a con- 
ventional (Euler) method vs. the single-plane method, 
for the specific case of rotations around the y axis. For 
the Euler method, if we have two fixed axes of rotation 
(x and z), the result corresponds to a three-step pro- 
cedure given by 7£(y, 0) = 7£(z, 7r/2)7£(x, 0)7£(z, — tt/2). 
In contrast, the single-plane method is accomplished in 
just two steps, with 7£(y, (j)) = 7£(n2, 7r)7£(z, 7r) where 
n2=(sin |, 0, cos |). Figure mc) shows an explicit con- 
struction of a 37r/2 rotation about the y axis, employing 
two rotations around axes in the xz plane. 



IV. ROTATIONS FOR EXCHANGE-ONLY 
QUBITS 

We now consider a concrete physical example: a logi- 
cal exchange-only qubit encoded in a three-spin block, for 
which all qubit operations are implemented using electri- 
cally tunable exchange couplings between the constituent 
spins, without requiring a magnetic field [6]. The effec- 
tive Hamiltonian for the spin system is given by 



where the exchange coupling parameters J12, J23, and 
J31 are typically non-negative. The total spin 5 t ot and 
its z component S£ ot are good quantum numbers since 
they commute with the Hamiltonian; we will use them 
to label the energy eigenstates. We are specifically inter- 
ested in the states with S to t= 1/2 and S* ot =±l/2, which 
are both two-fold degenerate. We specify these states as 
{ I £tot •> ^tot 5 } •> adopting the label /=0,1 for the degen- 
erate states. With these definitions, we can encode the 
qubit in a decoherence free subsystem as follows 0: 



|0)l 
|1)l 



1 1 

2'2' 
1 1 

2'2' 



1 1 

2'~2' 
1 1 

2'~2' 



(18) 
(19) 



where ||,|;0)=|5)i2®| t>3, \h \\ 1)=- V ^T/3\T ) 12 ® 
I t>3+x/273|T+) 12 ® I |) 3 , and ||,-|;0> and ||,-|;1> 
are defined as spin-flipped states of ||, |; 0) and ||, |; 1), 
respectively. 

In the logical qubit space, the Hamiltonian becomes 



H 



J12 + J23 + J31 „ . V3 ( J23 - ^31) 
4 1 + 4 a * 

— 2J12 + J23 + J31 



(20) 



The unitary operator U (t)=exip(—iHt/H) rotates the log- 
ical qubit around an axis n in the xz plane by an angle 
(/>, given by 

n = (VS ( J 23 - J31) , 0, -2 J 12 + J 23 + J31) , (21) 
(j) = Jt/h , (22) 



where J=^Jf 2 + ^23 + ^31 ^12^23 - J23J31 ^31^12- 
The decoherence free subspace considered in Ref. [6] cor- 
responds to the special case of a=l and b=0 in Eqs. (18) 
and (pi. 



Equation (17) suggests a ring-like coupling configura- 



tion for the quantum dots. Such configurations have been 
achieved in the laboratory [16l [17] ; however full control 
of the couplings is challenging. When full control is pos- 
sible, the qubit can be rotated along any axis in the xz 
plane. Single-qubit operations then follow the procedures 
described above. 

A more common experimental arrangement is the lin- 
ear triple quantum dot geometry, with one electron per 
dot. In this configuration, great tunability can be at- 
tained in the exchange couplings [18-23 . However, since 



H — J12S1 • S2 + J23S2 • S3 + J31S3 • Si 



(17) 



one of the exchange couplings in Eq. (17) is assumed to 
vanish, it is not possible to implement arbitrary rota- 
tions in the xz plane; only 2/3 of the plane is covered. 
For example, if J31 = 0, rotations are limited to the 
range tt/3 < < tt and < 6 < 2tt/3. Figure [3] shows 
the viable rotation axes in the xz plane when one of the 
exchange couplings is set to zero. 



(a) 



(b) 



(c) 



FIG. 3. (Color online) Regions of the xz plane where viable 
rotation axes are found, when one of the three exchange cou- 
plings is set to zero: (a) J31 = 0; (b) J23 = 0; (c) J12 = 0. 
Solid blue regions represent the rotation axes n, consistent 
with Eq. (21), assuming non-negative exchange couplings. 
Striped, red regions correspond to rotations around the — n 
axis by angle <fi, which is equivalent to rotations around n by 
angle 2ty — (j). 



Despite the fact that we cannot perform rotations over 
the entire xz plane for the linear dot geometry, many 
gates of interest can still be implemented in two or fewer 
steps. Focusing on the configuration J13 = 0, we see that 
rotations about the x axis (J23 = 2 J 12) an d the z axis 
( J23 = 0) can be accomplished in a single step. Rotations 
around the y axis by an angle (j) can be accomplished in 
two steps, by using one of the two solutions 

K(y, </>) = K(z, 7r)ft(ni, tt) = -ft(n 2 , tt)7J(z, tt) , (23) 

where ni=(sin f , 0, — cos ^) and n2 = (sin ^, 0, cos ^). Up 
to a global phase, the phase gate 5=[(1, 0), (0, i)] and the 
7r/8 gate T=[(l, 0), (0, exp(z7r/4))] can be implemented 
with single-step rotations around the z axis. 

The Hadamard gate if=l/>/2[(l, 1), (1, -1)] corre- 
sponds to a rotation around an axis with = 7r/4 and 
ip = 0, by an angle = tt. For the configuration with 
J12 = 0, this can be implemented in a single step. How- 
ever, in the J23 = or J31 =0 configurations, it can- 
not be implemented in fewer than three steps. If many 
Hadamard gates are required for a given quantum circuit, 
this could present a bottleneck. In this case, it might be 
preferable to change the logical qubit basis in Eqs. (18) 



and (19), to one where qubits 2 and 3 are exchanged. 



In the latter configuration, the Hadamard gate is accom- 
plished in one step. 

We can compare our single-plane rotation method with 
the serial gating scheme for exchange-only qubits, which 
was described in [6]. There it was shown that when only 
one exchange coupling (J12, «/23> or J31) is allowed at a 
time, then a general logical qubit rotation requires three 
(four) steps for a ring (linear) geometry, as illustrated in 
Figs. [4^ a) and (b). The single-plane rotation method de- 
scribe above can be viewed as a parallel gating scheme, 
where multiple exchange couplings are implemented si- 
multaneously, as illustrated in Figs. |4^c) and (d). As 
we have shown, when two (three) exchange couplings are 
allowed simultaneously, then a logical qubit rotation re- 
quires three (two) steps. Comparison with Figs. [4^ a) and 



(b) 



(d) 



FIG. 4. Logical qubit rotation schemes for exchange-only 
logical qubits. The horizonal lines depict the three physical 
spins comprising the logical qubit, and the vertical arrows 
indicate exchange couplings between two spins. The gating 
scheme of Ref. [6] considers serial exchange couplings in (a) 
a ring geometry, or (b) a linear geometry. The single-plane 
rotation scheme described here considers simultaneous, par- 
allel exchange couplings in (c) a ring geometry, or (d) a linear 
geometry. In the ring geometry, logical qubit gates require 
up to three steps in serial mode [(a)], and up to two steps in 
parallel mode [(c)]. In the linear geometry, logical qubit gates 
require up to four steps in serial mode [(b)], and up to three 
steps in parallel mode [(d)]. 



(b) shows that the single-plane method is always more ef- 
fecient than the serial gating scheme. 



V. CONCLUSIONS 

We have shown that the ability to tune qubit rotation 
axes in a single, fixed plane enables efficient, two-step 
implementations of single-qubit gates. This should be 
contrasted with fixed-axis methods (e.g., Euler angles), 
that require up to three steps. Our results can be ap- 
plied directly to logical qubits in the decoherence free 
subspaces and subsystems of three-spin logical qubits, 
such as exchange-only qubits. 

Our scheme can also be adapted to any qubit imple- 
mentation with at least partial control over the rotation 
axes. For example, the effective field corresponding to 
rotations of a singlet-triplet logical qubit [5] consists of a 
fixed B x component and a tunable, positive B z compo- 
nent. The resulting effective rotation axis covers about 
half of the xz plane [24 , 25 , thus enabling efficient gating 
methods similar to those described here. 

This work was supported by the DARPA QuEST pro- 
gram through a grant from AFOSR, by NSA/LPS 
through grants from ARO (W911NF-08-1-0482, 
W911NF-12-1-0607, and W911NF-09-1-0393), and 
by the DOD. The views and conclusions contained in 
this document are those of the authors and should not 
be interpreted as representing the official policies, either 
expressly or implied, of the US Government. 
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Appendix A: Rotation operators 

A rotation operator 7£(n, <p) is defined by its rotation 
axis n and rotation angle <fi. The rotation axis is de- 
scribed by a polar angle in the range [0, tt] and an az- 
imuthal angle i/j in the range [0, 2tt): 



n = (sin # cos?/;, sin # sin?/;, cos#) . 



(Al) 



The rotation angle is also in the range [0, 2tt). In terms 
of Pauli operators, the rotation operator can be written 

as 



exp 
1 cos 



icr • n sin - 



(A2) 



Any rotation operators in a two-dimensional Hilbert 
space can be represented as 



1Z = aot + i(aicr x + a 2 (7 y + a 3 cr 2 ) 
where ao,i,2,3 are all real. In matrix form, 



K = 



ao + ias a 2 + iai 
— a 2 + ia\ a — m 3 



(A3) 



(A4) 



The unimodular condition requires that ao+a 2 +a 2 +a 2 = 
1. 

The necessary and sufficient condition for two rotation 
operators 1Z\ = a^t + i{a\<j x + a^Oy + ^3^) an d 7£ 2 = 
b^ + i(bi<j x + b2<Jy + b3<j z ) to be identical, is that a$ = b$ 
for z = 0, 1, 2, 3. This is obvious by comparing the matrix 
elements of the two rotation operators. 



Appendix B: Two-step implementations of single 
qubit gates 

We want to identify two sequential rotation steps 
around axes in the xz plane that produce an arbitrary 
single-qubit gate, up to a global phase: 

ft(n, 0) = e^(n 2 , fcM&ufa). (Bl) 

According to the constraints discussed in the main text, 
n can be any direction, while hi and n 2 must lie in the 
xz plane: 

n = (sin # cos?/;, sin # sin?/;, cos#) , (B2) 
hi = (sin#i, 0, cos#i) , (B3) 
n 2 = (sin(9 2 ,0,cos6> 2 ) . (B4) 

We now transform the problem to a new set of (primed) 
coordinate axes, for which n lies in the new y' z' plane: 



n = sin 0'f + cos 0'z' 



(B5) 



This can be achieved by rotating the original coordinate 
axes around y axis. (See Fig. 2 in the main text.) Hence, 



the y axis is unaffected, while the x and z axes become 
x 7 and z' ' . The rotation angle (3 for this transformation 
is given by 



cos (3 = 



sin (3 = 



(B6) 



In the primed coordinate system, we easily find that ?// = 
7r/2 when n y > 0, and ?// = 37r/2 when n y < 0. 0' is 
obtained from 



sin 0' = sin sin ip . 



(B7) 



Expanding both sides of Eq. (Bl ) in terms of Pauli oper 



ators and matching their coefficients, we obtain the fol- 
lowing relations in the primed coordinate system: 



cos ■ 



VI //)/ n f \ • • VI 

cos — cos cos(c7o — 6-1 ) sin — sin — 

2 2 v 2 17 2 2 



= e 
k sin 0' sin 



sin #1 cos — sin h sin On sin — cos — 

1 2 2 2 2 2 



-sin(^ 2 -^)sin^sin^ 



cos c/ sin — = e 
2 



cos 6/1 cos — sin — 
1 2 2 



/)/ • 92 91 
- cos 6/9 sin — cos — 
2 2 2 



(B8) 
, (B9) 

(B10) 
(Bll) 



Here, we define /c = sgn[n y ] = sgn[sin#sin?/;]. Nu- 
merical results suggest that there will be a continuum of 
solutions for these equations. Here we set (/> 2 = tt, to sim- 
plify the equations and to enable an analytical solution. 
In this case, the equations reduce to 



cos I = -e^ cos(6> 2 - 6>0 sin ^ 
O = e^sin0'cos^ , 
k sin 0' sin ^ = -e ir? sin(<9 2 - 6>; ) sin ^ 



cos V sin - 



e i7? cos6> 2 cos^ 
2 2 



(B12) 
(B13) 
(B14) 
(B15) 



From Eq. (B13L we see that either sin On = 0, or cos ?r 
0. To satisfy Eq. (B15), we must have sin# 2 
therefore obtain 0' 2 = 0. 

The three remaining equations are 

\ 



k sin 0' sin 



-e ' cosc/i sin 



e ir > sin (9; sin ^ 
1 2 



cos c/ sm - 



0. We 

(B16) 
(BIT) 
(B18) 



From Eq. (B17), we see that e %v = k, since the sine func- 
tions are all positive for the range of angles 0', 0[ G [0, tt] 
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and 



e [0, 27r). Hence, 

/c cos — = — cos 6^ sin — , 



sintf sin — = sintfi sin 
2 1 



2 ' 



k cos ^ sin - = cos — 
2 2 



(B19) 
(B20) 
(B21) 



Here, we have two unknowns (0[ and 0i) and three equa- 
tions. However, the three equations are n ot indepen- 
dent. By squaring both sides of Eqs. ( |B19 ) and (B20) 
and adding them, we obtain 



cos — + sin 6/ sin — = sm — , (B22) 



which leads to 



cos sm — = cos — . 

2 2 



(B23) 



This is the same as the square of Eq. ( B21 ). We 
can obtai n another equation by dividing Eq. (B20) by 
Eq. (|B19l): 



k sin 0' tan | = - tan 0' x . (B24) 



(right-hand-side) of Eq. ( B19) as L\ (Ri), and similarly 
for L 2 (R 2 ) in Eq. dB2Q[ ), then Eq ( |B21| ) implies that 
+ ^2 = #i + and Ec l- d B24 l > leads to L 2 jL x = 
R 2 /Ri = 7. Here we note that L2 and i ?2 are both 
positive, as was explained below Eq. ( |B17 ). Plugging 
L\ L2/7 and R 1 R 2 /j into L\ + L% = R\ 
we obtain (I/7 2 + \)L\ = (I/7 2 + and then L 2 , 



^2? 



2 — 

Since L 2 and i?2 are positive, we obtain L 2 = R 2 . 
Now, from L 2 jL\ = R 2 jR\, we obt ain L i = R\. We 
can therefore determine from Eq. ( B21|) an d 0[ from 
Eq. ( |B24D . Note that Eqs. flB21| ) and QB24| ) uniquely 
determine 0i in the range [0, 2tt) and 6[ in the range 
[0,4 

To summarize, we can always implement an arbitrary 
single-qubit gate with two rotation steps around axes in 
the xz plane, gi ven by <p 2 = tt and 2 = 0, with <£ q ob- 



tained from Eq. ( |B21| ), and 0[ obtained from Eq. ( |B24 ). 
Of course, this is not the only solution. For example, we 
can also find a solution by choosing <pi = tt. In that case 
6[ =0, e %r] — — fc, and after similar procedure we obtain 



/c cos sm — = cos — , 
2 2 ' 



/c sin 7 tan — = tan 0. 



2 5 



(B25) 
(B26) 



We now s how t hat o nce \ and <pi are obtained, by 
solving Eqs. ( B21|) an d (B24), the results will also satisfy 
Eqs. (B19) and (B20). If we represent the left-hand-side 



which determine 2 and <p 2 uniquely. Once we determine 
6[ and 2 , we can transform back to the original coordi- 
nate system using 0i=0[ + j3 and 2 =0' 2 + /3. 
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